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PRECEDING  PAGE  BI^jWUNOT  EILMSD 


1.  INTRODUCTION  TO  FLAIR 


The  FLAIR  code  (acronym  for  FoLd  AIR)  was  developed  to  easily 
and  inexpensively  generate  x-ray  and  gamma-ray  transport  results  in  in- 
finite air  by  utilizing  an  existing  data  base  of  x-  and  gamma -ray  air  transport 
computed  by  Monte  Carlo.  FLAIR  provides  energy,  angle,  and  time  dependent 
photon  transport  results  from  energy  and  time  dependent  photon  sources  in 
infinite  air.  Results  are  possible  out  to  ranges  of  fifteen  mean -free -paths 
and  are  scaled  with  air  density. 

The  Monte  Carlo  calculations  were  performed  with  source  energy 
bands  so  that  arbitrary  source  spectrums  can  be  used.  Time  convolution 
is  also  performed  for  time  dependent  sources,  although  the  energy  and  time 
dependence  of  a source  must  be  separable.  The  photon  transport  data  base 
is  stored  on  two  magnetic  tapes,  one  for  x-rays  and  one  for  gamma -rays. 

The  x-ray  tape  has  22  energy  bands  from  10  keV  to  300  keV  and  the  gamma- 
ray  tape  has  25  energy  bands  from  10  keV  to  14  MeV.  Calculations  below 
300  keV  can  be  performed  with  the  gamma -ray  tape,  although  better  energy 
resolution  would  be  obtained  with  the  x-ray  tape. 

Problem  specification  is  very  flexible  allowing  the  user  to  choose 

. 

from  a variety  of  units  and  input  and  output  options.  For  example,  a user 
could  specify  a source  in  photons  AeV  and  generate  forward  energy  current 
in  calories  /keV,  or  the  source  could  be  MeV/group  and  the  results  in 
rad(Si)/sec.  This  user  flexibility  will  alleviate  most  of  the  pre-processing 
and  post -processing  of  data  associated  with  a FLAIR  run. 
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The  x-ray  and  gamma -ray  transport  data  base  was  generated  by 
the  authors  with  the  DART  Monte  Carlo  code  (documented  in  Ref.  1)  which 
was  developed  specifically  for  this  task.  The  DART  code  uses  point  cross 
sections  and  rigorously  samples  from  the  coherent  and  incoherent  scattering 
angular  distributions.  Flux  or  current  quantities  computed  with  DART  can 
be  a joint  function  of  up  to  five  independent  variables:  two  spatial,  energy, 
time,  and  angle.  Although  DART  generates  a range  and  angle  position  de- 
pendence from  a "gun  barrel"  source,  FLAIR  treats  only  isotropic  sources, 
having  integrated  out  the  position  angle  dependence.  A report  documenting 
the  data  base  and  discussing  our  investigation  of  coherent  and  incoherent 
scattering  effects  will  be  published. 

Since  single  Monte  Carlo  calculations  were  used  to  generate  data  at 
all  ranges,  energies,  times  and  angles,  our  general  guideline  for  statistical 
accuracy  was  to  have  less  than  10%  fractional  standard  deviations  for  the 
energy  spectrum,  time  dependence,  and  angle  dependence  at  all  ranges.  This 
was  most  difficult  at  the  deep  penetrations  and  the  late  times  (especially  at 
close-in  radii).  Importance  functions  and  running  times  were  set  so  that 
this  10%  accuracy  was  met  in  most  instances. 

Several  important  features  distinguish  the  FLAIR  code  from  the  ATR 
code;  these  are  summarized  in  Table  1.  First  ATR  is  more  general  than 
FLAIR  including  neutrons  and  secondary  gamma  rays  as  well  as  prompt  x- 
rays  and  gamma -rays.  FLAIR  however  provides  much  more  comprehensive 
data  for  photons  that  includes  full  energy,  time,  angle  and  range  dependence. 
FLAIR  results  are  also  more  accurate:  first,  the  the  data  base  calculations 
with  DART  used  point  cross  sections  instead  of  the  multigroup  approximation 
used  by  ATR;  and  secondly,  FLAIR  uses  the  raw  data  without  the  intermediate 
step  of  curve  fitting.  Notice  that  both  codes  use  the  free-field  input  and  a 
variety  of  output  options.  ATR  also  allows  air -over -ground  and  high  altitude 
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TABLE  1.  COMPARISON  OF  THE  FLAIR  AND  ATR  CODES 


Property 

FLAIR 

ATR 

Data  Base 

point  cross  section 
Monte  Carlo  with  C-I 
scattering,  x-ray  and 
gamma -ray  only 

multigroup  discrete 
ordinates  and  Monte 
Carlo,  K-N  scattering, 
also  neutron  and 
secondary  gamma-ray 

Data  Acquisition 

directly  from  data  base 
tapes 

from  parameterizations 
(effects)  of  data  base 

Data  Content 

: i 

complete  range,  energy, 
time,  angle  dependent 
flux 

no  time  dependence  for 
x-rays,  range,  energy, 
time  flux  or  range, 
energy,  angle  fluence 
for  gamma -rays 

Convenience 

Features 

free -field  input,  variety 
of  output  options 

free -fie  Id  input,  variety 
of  output  options,  air- 
over -ground  and  high 
altitude  correction  factors, 
integral  mass  scaling, 
built-in  sources  and 
responses 

Approximate 
time /run 

15  seconds 
CDC  7600 

15  seconds 
UNI  VAC  1108 

Approximate 

cost/run 

$5.00 

$.  50 

correction  factors  which  are  not  present  in  FLAIR.  For  these  reasons  the 
FLAIR  code  can  provide  a far  more  complete  and  accurate  picture  of  photon 
transport  in  a one -dimensional  atmosphere  than  ATR.  On  the  other  hand, 
with  commensurate  sacrifices  in  accuracy,  ATR  is  faster  running  and  pro- 
vides transport  data  for  all  burst  and  target  altitudes.  ATR  is  oriented 
toward  systems  and  "quick  look"  studies  while  FLAIR  is  oriented  toward 
benchmarking  radiation  environments  from  weapon  systems  and  calculating 
differential  fluxes  for  input  to  radiation  response  codes. 
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2.  DESCRIPTION  OF  FLAIR 


This  chapter  describes  in  detail  the  workings  of  FLAIR.  The  format 
of  the  data  base  is  presented  along  with  equations  for  all  operations  per- 
formed by  FLAIR.  As  each  function  of  FLAIR  is  described,  the  correspond- 
ing input  card  formats  are  presented  in  detail.  For  user  convenience  these 
input  instructions  are  summarized  in  Appendix  A. 

2.  1 FORMAT  OF  THE  DATA  BASE 

The  x-ray  and  gamma -ray  tapes  contain  the  radiation  transport  data 
in  the  form  of  4ffR  photon  currents  per  group  as  a function  of  range,  angle, 
time,  energy  and  source  energy.  This  data  has  not  been  parameterized,  and 
the  tapes  contain  one  number  for  each  five  dimensional  value.  The  format  of 
the  tapes  is  presented  in  Table  2.  Notice  that  the  order  of  the  variables 
from  the  fastest  to  the  slowest  is  energy,  time,  angle,  range  and  source 
energy,  and  that  the  zeros  for  energies  higher  than  the  source  energy  are 
not  written  on  the  tape.  Table  3 presents  the  energy  bounds  of  the  22  x- 
ray  and  25  gamma -ray  energy  groups.  Because  the  gamma -ray  tape  covers 
the  x-ray  energies  as  well,  it  can  also  be  used  for  x-ray  calculations.  How- 
ever more  accurate  energy  resolution  will  be  obtained  by  using  the  x-ray 
tape  for  x-ray  problems.  These  energy  bands  presented  in  Table  3 were 
used  for  both  the  source  energy  bands  and  the  detector  energy  scoring  bins. 

The  seventeen  ranges  for  both  x-rays  and  gamma -rays  included  in 
the  data  base  were  chosen  as  a number  of  mean -free -paths  of  the  source 
energy.  The  data  is  presented  as  currents  exactly  at  these  ranges  as  cal- 
culated with  a boundary  crossing  estimator.  The  first  two  radii  are  at 
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TABLE  2.  FORMAT  OF  THE  FLAIR  DATA  TAPE 


Record  1 - 224  words 

Hollerith  header  record  giving  references  of  the  documents  de- 
scribing the  work. 

Record  2 - 600  words 


Header  record  for  each  source  energy  band. 

IS0R  - source  energy  number 

NR  - number  of  radii 

NE  - number  of  energy  bins 


- number  of  time  bins 


NMU  - number  of  angle  bins 

LENGTH-  number  of  data  for  this  source  band 
= NR*(NE  -IS0R+1)  *NT*NMU 

LBUFF  - record  length  of  each  tape  record  for  the  data  = 5000 

NREC  - number  of  records  necessary  to  transmit  LENGTH 
words.  The  last  record  may  be  a partial  record. 

ES0RHI  - upper  energy  of  source  band  (rest  mass  units) 

ES0RL0  - lower  energy  of  source  band  (rest  mass  units) 

RBIN(18)-  scoring  radii  (these  are  different  for  each  source  band 

as  each  DART  calculation  went  out  to  15  mean  free  paths) 

EBIN(25)-  energy  bin  boundaries  (rest  mass  units) 

TBIN(18)-  time  bin  boundaries  (seconds) 

UBIN(6)  - cosine  bin  boundaries 

UNC0LD(18)  - uncollided  contribution  at  each  radius 

SIR(18)  - fractional  standard  deviation  of  the  total  current  at 
each  radius 


NREC 


TABLE  2.  (Cont'd.) 


SIRE(450)  - fractional  standard  deviations  of  the  energy  spectra 

total  current  at  each  radius.  Energy  is  the  fastest 
varying  index. 

Record  3 -5000  words  (or  LENGTH  if  NREC  = 1)  of  four  dimensional  dif- 
ferential currents.  The  order  of  variation  of  indices  from  fastest 
to  slowest  is  NE,  NT,  NMU,  NR. 

Record  3 is  repeated  until  the  entire  source  band  is  completed. 

Records  2 and  3 are  then  repeated  until  all  source  bands  are 
completed. 


n 


TABLE  3.  FLAIR  X-  AND  GAMMA -RAY  ENERGY  BINS 


Bin 

gamma -ray  upper 
energy  bounds 

x-ray  upper 
energy  bounds 

1 

14. 0 MeV 

300.  keV 

2 

10.0 

260. 

3 

8.0 

220. 

4 

7.0 

190. 

5 

6.0 

160. 

6 

5.0 

140. 

7 

4.0 

120. 

8 

3.5 

105. 

9 

3.  0 

90. 

10 

2.5 

80. 

11 

2.0 

70. 

12 

1.66 

65. 

13 

1.33 

60. 

14 

1.0 

55. 

15 

800.0  keV 

50. 

16 

600.0 

45. 

17 

450.0 

40. 

18 

300.0 

35. 

19 

220.0 

30. 

20 

160.0 

25. 

21 

120.0 

20. 

22 

90.0 

15. 

23 

65.0 

10.  lower  bound 

24 

45.0 

25 

30.0 

10. 0 lower  bound 
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. 25  and  .66  mean-free -paths  and  the  next  15  radii  go  from  1 to  15  mean-free- 
paths,  respectively.  Because  the  ranges  depend  on  the  source  cross  section, 
a calculation  with  a source  energy  spectrum  will  require  a different  range 
index  and  interpolation  value  for  each  source  energy  in  the  spectrum.  Table 
4 presents  the  17  time  bins  and  six  angle  bins  used  in  both  the  x-  and  gamma- 
ray  tapes. 

The  short  initial  record  contains  the  boundaries,  the  uncoil ided  contri- 
butions and  fractional  standard  deviations  of  a few  quantities  in  the  data  base. 
For  angle  folding  and  time  convolution,  FLAIR  treats  the  uncollided  flux  as  a 
seventh  angle  bin  and  an  18th  time  bin.  Thus  when  converting  current  to  flux, 
the  uncollided  can  be  given  its  correct  response  of  1.0  rather  than  the  most 
forward  angle  bin’s  value  of  1.2.  Likewise  when  performing  the  time  convolution, 
the  uncoil ided  is  treated  as  a delta  function  contribution  at  time  zero  rather 
than  averaged  over  the  first  time  bin.  At  present  FLAIR  makes  no  use  of  the 
statistical  information  written  on  the  first  record. 

2.  2 GENERAL  FORM  OF  INPUT  CARDS 

(31 

The  input  routines  for  FLAIR  were  borrowed  from  the  ATR  code' 
and  as  such  have  the  same  general  format.  Problem  specifications  are 
input  in  free -field  format  on  various  types  of  cards.  An  asterisk  (*)  in- 
dicates the  start  of  a new  card  type.  A character  or  characters  following 
the  asterisk  indicate  which  type  data  will  follow.  For  example  ranges  are 
input  on  *R  cards,  air  density  on  a *DEN  card,  and  the  energy  spectrum 
on  a *ESOR  card.  Usually  the  units  are  specified  after  the  type.  For 

ranges,  units  on  the  *P.  v-id  may  be  KM,  KFT,  MI,  FT,  M,  CM,  or  GM 
2 

(for  grams/cm  of  air).  After  the  units,  the  values  are  specified  on  the 
remainder  of  the  card  or  continued  on  further  cards  without  an  asterisk. 

Some  card  types  do  not  contain  values,  but  merely  specify  output  or  ex- 
ecution options.  All  input  cards  will  be  explained  in  the  following  sections 
which  describe  the  FLAIR  operations. 
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TABLE  4.  FLAIR  TIME  AND  ANGLE  BINS 


Bin 

time  bin 

angle  bin 

lower  bound 

upper  cosine 

1 

0.0 

1.0 

2 

7. 0(-9)  sec 

.666 

3 

1.5(-8) 

.333 

4 

3. 0(-8) 

0.0 

5 

6. 0(-8) 

-.333 

6 

1*  2(— 7) 

-.666 

7 

2.  4( -7) 

-1.0  lower  bound 

8 

3.  5( — 7) 

9 

6. 0( — 7) 

10 

9.0(-7) 

11 

1.2(-6) 

12 

1.6(-6) 

13 

2. 0(-6) 

14 

2.  5(-6) 

15 

3.  2(-6) 

16 

5. 0(-6) 

17 

9. 0(-6) 

1.  5(-5)  maximum  time 
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The  list  of  values  on  a FLAIR  (or  an  ATR)  card  is  used  to  specify 
the  numerical  data  a card  may  require.  A number  in  a list  may  appear  in 
a variety  of  forms  to  suit  the  particular  user  or  problem.  For  example, 
some  of  the  forms  in  which  the  number  400.  may  appear  are  as  follows: 

400  400  . 4.E+2  4E+2  4+2  4000-1 

At  least  one  number  must  appear  in  a list  of  values  . Two  or  more 
numbers  are  separated  from  one  another  by  the  occurrence  of  one  or  more 
blank  characters.  Therefore,  the  user  is  restricted  from  specifying  a 
number  in  which  internal  blanks  appear,  or  which  is  split  on  two  or  more 
card  images.  A further  restriction  exists  upon  the  magnitude  of  such 
numbers:  since  a number  is  interpreted  as  a function  of  up  to  three  integer 
parts  (a  whole  part,  a fractional  part,  and  an  exponent),  none  of  the  parts 
of  a number  may  exceed  in  magnitude  the  greatest  integer  value  appropriate 
to  the  host  machine,  nor  can  the  number  generated  from  these  three  parts 
exceed  the  host  machine’s  allowable  range  of  representable  numbers. 

2.  3 RANGE  INTERPOLATION  AND  DENSITY  SCALING 

Values  can  be  calculated  by  a single  FLAIR  eAecution  for  up  to  20 
ranges.  These  ranges  are  input  on  the  *R  card  which  has  the  format 

*R  units  values 
where  the  units  may  be 


KM 

- kilometers 

KFT 

- 1000  feet 

MI 

- miles 

FT 

- feet 

15 
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M - meters 

CM  - centimeter  (default  option  used  if  units  left  blank) 

GM  - gram/cm^  of  air 

More  than  one  *R  card  may  be  used  with  different  units,  but  the 
maximum  number  of  ranges  per  run  is  20.  At  least  one  *R  card  is  required 
for  each  execution. 

When  the  GM  units  option  is  used,  the  density  known  by  FLAIR  when 
the  *R  card  is  read  is  used  to  set  the  distance  in  centimeters.  Therefore 
the  density  should  not  be  changed  by  having  a *DEN  card  after  the  *R  card. 

If  a density  other  than  1. 29  x 10  "3  gm/cm3  (sea  level)  is  desired,  it  should 
be  set  by  a *DEN  card  before  any  *R  GM  cards  are  read. 

Values  are  then  obtained  from  the  data  base  by  semi -logarithmic 
interpolation  on  the  ranges.  If  results  are  desired  at  RCALC,  the  data 
base  is  searched  for  the  values  at  the  radii  immediately  less  and  greater 
than  RCALC.  These  values  are  stored  respectively  in  two  partitions  of 
blank  common  (referred  to  as  BClegg  and  BC  and  later  interpolated 

into  a third  partition,  BCgum.  The  semi -logarithmic  interpolation  equations 
then  are 


where 

, 

I . 


If  RCALC  is  greater  than  the  maximum  radius  in  the  data  base  (15  mean -free - 
paths),  then  BCsum  = 0 is  assumed.  Likewise  if  RCALC  is  less  than  the  mini- 
mum radius  (.25  mean-free -paths),  then  BCgum  = BC  tej>  is  assumed  for 


BC  = BCa  . Be!1"®) 

sum  greater  less 


Rgreater  " RCALC 
a-  « 


greater 


less 
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all  values  except  the  uncollided  which  is  interpolated  between  1.0  and  the  un- 
collided at  the  first  radius. 

FLAIR  works  by  reading  the  data  base  tape  and  storing  or  summing 
the  data  in  either  a less  than  (BClegg)  or  a greater  than  (BCgreater)  array 
in  blank  common  when  appropriate.  After  all  the  data  for  a single  source 
energy  have  been  processed,  the  two  arrays  are  interpolated  and  summed 
into  a third  array  (BCgum)  which  contains  the  running  sum  over  all  source 
groups  after  interpolation.  BClegg  and  BCgreater  are  zeroed  before  read- 
ing the  next  source  band. 

FLAIR  results  can  be  scaled  by  air  density  by  including  a density  on 
a *DEN  card  whose  format  is 

*DEN  units  value 

where  the  units  may  be 

O 

GM/CM3  - grams/cm  (default  if  no  units  are  specified) 
AT/BRNCM  - atoms /barn-cm 

O 

If  no  density  card  is  encountered,  the  sea  level  density  of  1.  29  x 10 
3 

gm/cm  is  assumed.  The  data  base  calculations  were  performed  at 
-3  3 

oQ  = 1.00  x 10  gm/cm  air  density.  To  scale  the  FLAIR  results  all 

input  radii  are  multiplied  by  (p/pQ)  to  generate  calculation  radii  (RCALC). 

The  RCALC  values  are  used  for  range  interpolation  while  the  input  radii 

2 

are  used  to  divide  the  final  results  by  4^R  if  required.  The  time  dependence 

is  also  scaled  by  multiplying  the  time  bin  boundaries  by  (p  /p).  In  this 

2 ^ 

manner  values  like  photon/cm  -sec  will  be  divided  by  (p/pQ)  , two  factors 
from  the  radii  and  one  from  the  time  bins. 

Unlike  the  ATR  code,  FLAIR  does  not  have  the  air  density  table 
built  in  and  as  such  cannot  perform  integral  mass  scaling  automatically. 

However  this  can  be  done  by  hand  by  inputing  the  ranges  with  the  GM 

2 2 
(grams/cm  ) units  option  on  the  *R  card  and  requesting  4ttR  results 
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on  the  *EXC  card  (to  be  explained  later).  These  results  can  then  be  divided 
2 

by  the  real  4trR  giving  correct  results  for  sources  and  detectors  at  dif- 
ferent altitudes. 


2. 4 FOLDING  AND  OUTPUT  OPTIONS 

FLAIR  allows  the  results  to  be  folded  with  an  arbitrary  energy  re- 
sponse and  an  arbitrary  angular  response.  Two  energy  responses  are  built 
in;  photons  and  energy  (in  several  units).  The  photon  response  is  all  unity 
and  the  energy  response  is  the  average  energy  of  the  group  in  the  desired 
units.  Three  angle  responses  are  built  in;  they  are  current,  forward  cur- 
rent, and  flux.  Since  the  raw  data  is  a current,  the  current  response  is 
all  1.0,  the  forward  current  response  is  three  1.0's  and  three  zero's  and 
the  flux  response  is  one  over  the  average  angle  in  an  angle  bin  which  is 
1.2,2.,  6. , 6. , 2. , and  1. 2.  An  angle  response  of  1. 0 is  used  to  convert 
the  uncollided  from  current  to  flux.  These  angle  and  energy  response  fold- 
ings are  performed  before  the  raw  data  is  stored  or  summed  in  BClegg  or 

BC  greater* 

The  user  can  also  input  his  own  energy  or  angle  response  with  either 
a *ERESP  or  a *ARESP  card.  The  format  for  the  *ERESP  card  is 

♦ERESPn  ♦title*  values 

n - 1 to  6 (maximum  of  6 energy  response  cards)  used  to 
identify  the  response. 

title  - a maximum  of  20  characters  to  be  printed  on  output. 

If  omitted  ERESPn  will  be  used. 

values  - one  value  needed  for  each  energy  bin  from  high  energy 
to  low  energy. 
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A *ERESPn  card  is  required  if  ERESPn  or  DERESPn  appears  on  a 
*0,  *T,  or  *P  card.  Similarly,  the  format  for  the  *ARESPn  card  is 


♦ARESPn  ♦title*  values 

- 1 to  7 (maximum  of  7 angle  response  cards)  used  to 
identify  the  response  in  the  code. 

- maximum  of  20  characters  used  to  identify  the  re- 
sponse on  the  output.  If  omitted  ARESPn  will  be 
assumed. 

- one  value  for  each  angle  bin  from  forward  to  back- 
ward, plus  one  more  for  the  uncollided.  If  the 
uncollided  value  is  not  included  it  will  be  set  to  the 
first  value. 

This  card  is  required  when  ARESPn  or  DARES Pn  is  specified  on  a *0,  *T, 
or  *P  card. 

FLAIR  allows  several  output  options  which  include  all  combinations 
of  energy,  angle,  and  time  dependence  at  the  desired  ranges.  In  other 
words,  the  results  may  be  integrated  over  any,  all,  or  none  of  the  angle, 
energy  and  time  variables.  Also  final  results  may  be  per  group,  normalized 
per  group,  differential,  or  normalized  and  differential.  All  these  options, 
as  well  as  the  desired  response  functions  are  specified  on  a *0  card.  The 
format  of  the  *0  card  is  as  follows: 


n 

title 

values 


m 


OUTPUT  CARD 


♦O  angle  energy  time  type  flag 

These  parameters  may  occur  in  any  order. 


angle  may  be  or 

FLUX 

FCUR 

CUR  (default) 

ARESPn  (7  is  max.  value  of  n) 


energy  may  be  or 


PH  (default) 

MEV 

KEV 

CAL 

ERESPn  (6  is  max.  value  of  n) 


time  may  be  or 

TIME  (default) 

type  may  be 

G for  group  data  output 

GN  for  normalized  group  data  output 

D for  delta  data  output  (default) 

DN  for  normalized  delta  data 

flag  may  be  blank  or  ONLY. 


L 


angle  dependent 

DFLUX 
DFCUR 
DC  UR 
DARESPn 


energy  dependent 

DPH 

DMEV 

DKEV 

DCAL 

DERESPn 


time  dependent 

DTI  ME 
DTSOR 
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At  least  one  *0  card  (or  *P  or  *T  card)  is  required  for  each  problem 
execution.  Notice  that  the  angle,  energy  and  time  parameters  specify  both 
the  response  to  be  used  with  that  variable  and  whether  that  variable  is  to  be 
summed  over  or  output  with  dependence.  Thus  FLUX  causes  the  angle  to  be 
summed  over  with  the  flux  response  while  DFLUX  causes  the  angle  dependence 
to  be  output  as  a flux.  The  equations  then  are 

FLUX  and  DFLUX 

7 

O = z ARESP  (flux).  F.  = ARESP  (flux).  F. 

where  F.^  is  a quantity  on  the  data  tape.  O and  O then  are  values  gene- 
rated in  the  BC^egg  or  BC  a^er  arrays.  Similar  equations  for  energy  and 
time  also  exist.  For  example  if  ERESP1  is  a silicon  flux  to  dose  conversion 
(rad(Si)/(y/cm2)), 

*0  FLUX  E RES  Pi  DTI  ME 

generates  the  time  dependent  silicon  dose.  The  summation  performed  would 
be 

NA  NE 

O,  = I E ARESP  (flux).  ERESPl.  F... 
k i=l  j=l  1 J «k 

k = 1,  NT  . 

For  a second  example 

*0  DFCUR  DCAL 

generates  the  time  integrated  forward  current  energy  spectrum  in  units  of 
calories  with  the  angle  dependence.  Notice  that  the  default  value  of  TIME 


may  be  left  blank.  This  summing  operation  generates  a two  dimensional 
array  which  looks  like 

NT 

O.,  = E ARESP  (fcur),  ERESP  (cal).  F... 

*J  jj=i  1 J 

i = 1,  NA  , j = 1,  NE  . 


The  type  parameter  on  the  *0  card  may  be  used  to  modify  the  output 
units.  The  G option  causes  the  results  to  be  output  as  generated  in  FLAIR 
on  a per  group  basis  without  dividing  by  any  bin  widths.  The  GN  option 
causes  any  variable  dependence  to  be  normalized  to  1.0  so  that  the  sum 
over  the  entire  array  (it  may  be  1,  2,  or  3 dimensional)  will  be  1.0.  The 
D option  (which  is  the  default  value)  for  delta  data  causes  the  unnormalized 
values  to  be  divided  by  all  appropriate  bin  widths.  For  angle  the  division 
is  per  cosine  not  per  steradian.  Thus  the  total  integral  over  a D option  re- 
sults would  equal  a summation  over  the  same  G option  results.  The  last  option 
DN  causes  the  data  to  be  first  normalized  and  then  divided  by  all  bin  widths. 

Normally  the  code  also  computes  all  lower  level  operations  as  well, 
so  that  if  the  user  requests  a three-dimensional  array  he  will  also  get  three 
two-dimensional  arrays,  three  one -dimensional  arrays,  and  a total  value. 

Both  printout  and  core  storage  can  be  minimized  if  ONLY  is  included  on 
the  *0  card.  In  this  case  only  the  operation  requested  will  be  performed, 
and  the  lower  order  operations  will  be  suppressed. 

Two  other  cards,  the  *T  and  the  *P  cards,  are  identical  to  the  *0 
card  except  that  they  also  cause  the  results  to  be  written  to  tape  or  punched 
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respectively  as  well  as  printed.  See  the  description  of  subroutine  OUT2 
in  Section  3. 2 for  a description  of  the  tape  and  punch  formats. 

2.  5 SOURCE  SPECTRUM  AND  YIELD  NORMALIZATION 

The  source  energy  spectrum  is  specified  on  a *ESOR  card  which 
has  the  following  format 

*ESOR  units  values 

units  - PH  - photons  /group  (default) 

MEV  - MeV/group 
KEV  - keV/group 
CAL  - calories /group 
PH/MEV  - differential  photons /MeV 
PH/KEV  - differential  photons /keV 

E/E  - either  MeV/MeV  or  keVAeV 
CAL/KEV  - calories  A^V  differential 

Values  are  from  the  highest  energy  group  to  the  lowest  energy  group.  If 
fewer  are  given  than  required,  the  lower  energies  will  be  assumed  zero. 
One  *ESOR  card  is  required  for  each  problem  execution. 

The  multiplication  by  the  energy  source  is  performed  at  the  same 
time  as  the  range  interpolation.  After  all  data  for  a single  source  band  is 
stored  or  summed  in  the  EClegg  and  Bcgreater  arrays  (and  the  angle  and 
energy  responses  are  included),  the  data  is  interpolated  on  range,  multi- 
plied by  the  energy  source  intensity  and  summed  into  the  BC  array. 
r s um 

The  number  of  photons  can  be  specified  on  a yield  card  or  *Y.  The 
format  for  the  *Y  card  is  as  follows: 


*Y  units 

value 

units  - 

PH  - photons  (default) 

KT  - kilotons  of  photons 
MEV  - total  energy  in  MeV 
CAL  - total  energy  in  calories 
MT  - megatons  of  photons 
MO  - moles  of  photons 

value  - 

one  value. 

If  the  *Y  card  is  omitted  the  source  spectrum  is  determined  from 
the  normalization  on  the  *ESOR  card.  If  the  *Y  card  is  included,  the  nor- 
malization is  set  by  the  *Y  card  and  any  normalization  that  exists  on  the 
*ESOR  card  is  ignored.  The  number  of  photons,  the  average  energy  of 
the  photons  and  the  yield  of  the  photons  in  kilotons  is  always  printed. 

2.  6 TIME  CONVOLUTION 


Time  convolution  of  a time  dependent  source  can  be  performed  in 
FLAIR  if  the  energy  dependence  and  the  time  dependence  of  the  source  are 
assumed  separable.  Clearly  time  convolution  is  only  important  if  the  time 
dependence  is  desired.  Therefore  time  convolution  is  signaled  by  DTSOR 
on  the  *0  card  (or  *P  or  *T  cards).  The  data  for  the  time  convolution  is 
input  on  a *TSOR  card  which  has  the  format 


*TSOR  (tk,wk,k=l,n) 

where  t^  and  w^  are  time  and  weight  pairs,  t^  has  units  of  seconds  and 
w^  has  units  of  photons.  Figure  1 shows  how  w^  corresponds  to  the 
number  of  photons  that  can  be  assumed  at  the  exact  time  t^.  Therefore 
wk  is  equal  to  the  area  of  the  rectangle  which  is  centered  at  t^  so  that 
the  time  dependence  is  assumed  to  be  n delta  functions  with  corresponding 
weights.  The  code  normalizes  the  w^  so  that 
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1 


TIME  (SEC) 


Figure  1 

' 


r 


Representation  of  the  Time  Dependent  Source  for  Time  Convolution. 
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n 

£ 


k=l 


wk = l-° 


The  following  discussion  describes  how  the  time  convolution  works 


in  FLAIR.  Let  T.  be  the  present  time  bin  boundaries  which  go  from  0 to 


Tn  and  let  o.  be  the  time  dependent  value  in  the  time  bin  between  T.  1 


and  T..  The  o.'s  are  calculated  in  FLAIR  as  if  the  *0  card  had  been 


DTIME  instead  of  DTSOR.  From  the  T/s,  t^’s  and  w^'s  a convolution 


matrix  M. . can  be  defined  such  that 


N 

C.  = £ M 

j=l 


. . (D  • 

1J  J 


where  the  C.'s  are  the  desired  convoluted  time  results.  Notice  that  M 
1 iJ 

is  independent  of  o.  so  that  the  convolution  matrix  M..  can  be  used  on 

J iJ 

several  sets  of  Oj  's*  can  be  further  broken  into  a sum  of  simpler 


matrices  by 


Mij  - j wk  mij 


Ir 

where  the  matrices  m..  depend  only  on  the  T.'s  and  the  one  time  t,  . 
k . ..  1 J . k 


Thus  m!\  is  the  convolution  matrix  for  a delta  function  source  at  a time 
other  than  zero. 


For  simplicity  the  k superscript  is  now  dropped  so  that  m..  is 
defined  as  a function  of  the  T.  ’s  and  the  time  t. 


* 

Using  linear  interpolation, 


it  can  be  seen  that  m. . is  defined  by 

ij  * 


where  T.  are  the  same  time  bins  as  T\.  These  equations  can  be  seen  by 

examining  the  overlap  of  two  sets  of  time  bins  which  are  identical  except 

that  one  is  displaced  by  t.  These  expressions  define  the  im.'s  which  can 

be  summed  to  define  the  total  convolution  matrix  M...  The  time  convolution 

U 

operation  is  performed  in  FLAIR  after  the  non-convoluted  time  dependence 
has  been  calculated  and  after  all  source  energy  bands  have  been  completed. 


2.  7 EXECUTION  CARDS 

Two  additional  cards  need  to  be  mentioned:  the  *EXC  card  and  the 
**  card.  The  *EXC  card  or  *EXC  4PIRSQ  signals  that  the  input  is  fully 
specified  and  FLAIR  begins  reading  the  data  tape  and  calculating  results. 
The  4PIRSQ  option  causes  all  output  results  (except  normalized  values)  to 


r 


1 


be  multiplied  by  4ffR  . If  no  4PIRSQ  appears  on  the  *EXC  card  the  results 
2 

are  per  cm  . After  FLAIR  finishes  the  problem,  the  data  tape  is  rewound 
and  FLAIR  tries  to  read  more  input  to  specify  another  problem.  The  ** 
card  terminates  the  job  on  a normal  exit. 

2.  8 CORE  STORAGE 

As  implied  in  this  section  but  not  yet  stated  the  number  of  operations 

performed  in  a single  calculation  is  limited  by  computer  printout  and  blank 

common  length.  For  each  number  printed,  punched,  or  written  to  tape, 

three  storage  locations  are  required  in  blank  common  for  BCless> 

and  BC . In  addition  5000  locations  are  used  in  BC  to  read  the  data  tape 

sum 

records.  Presently  the  entire  code  fits  in  30,00010  locations  with  blank 
common  set  at  20, 000  This  would  allow  one  full  three-dimensional  op- 
eration with  all  lower  operations  as  well.  Blank  common  however  can  be 
reduced  or  increased  to  fit  FLAIR  on  smaller  or  larger  computers,  with  a 
corresponding  change  in  number  of  operations  which  could  be  performed. 


3.  DOCUMENTATION  OF  FLAIR 


This  chapter  documents  the  subroutines  and  variables  used  in  the 
FLAIR  code,  and  as  such  will  be  useful  to  programmers  that  are  using  or 
modifying  FLAIR.  The  first  section  describes  all  variables  used  in  the 
named  commons  and  the  second  section  documents  the  subroutines.  Several 
routines  used  in  the  free-field  input  were  borrowed  from  ATR  and  used  as 
is  by  FLAIR.  As  such  they  are  not  documented  in  this  report.  It  should  be 
noticed  that  the  description  of  subroutine  CARDIN  gives  additional  informa- 
tion about  the  input  including  discussion  of  errors  in  input. 

3.  1 DESCRIPTIONS  OF  FLAIR  COMMON  BLOCKS 


Blank  Common  - Input  and  output  data  arrays 


BUF(5000)  - Contains  one  record  of  four  dimensional  differential  flux  data 
from  FIAIRIT  tape.  The  order  of  variation  of  indices  from 
fastest  to  slowest  is  NEBIN,  NTBIN,  NUBIN,  NR  BIN. 


BC(  15000)  - This  array  is  divided  into  three  parts,  BC,egs,  BC^^^, 

and  BCsum.  The  length  and  composition  of  these  subdivisions 
vary  with  the  number  and  type  of  folding  operations  requested, 
as  well  as  with  the  number  of  RCALC  requested.  The  values 
associated  with  RBIN  immediately  less  than  RCALC  are  sum- 
med in  BCiess.  The  values  associated  with  the  RBIN  immedi- 
ately larger  than  RCALC  are  summed  in  BCgreater.  The 
interpolated  results  between  BCiess  and  BC greater  are  ac~ 
cumulated  in  BCsum.  Each  of  these  three  arrays  are  further 
subdivided  into  NR  arrays,  one  associated  with  each  RCALC. 
Each  of  these  is  further  subdivided  into  NIOP  arrays,  one  for 
each  folding  operation.  The  storage  requirements  and  layout 
for  each  type  of  fold  array  is  described  in  Table  5. 
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'rom  fastest  to  slowes 


/ACCSRY/  - Hollerith  constants  for  input  processing  subroutines. 

IDP  - decimal  point. 

IEX  - letter  E. 

MIN  - minus  sign. 

IPL  - plus  sign. 

LPAR  - left  parenthesis. 

IBLK  - blank. 

IRPAR  - right  parenthesis. 

ISTR  - asterisk. 

/BIGGER/  - Maximum  storage  limit  for  BC. 

LIMIT  - maximum  storage  size  for  BC.  By  changing  this  value 

and  changing  the  blank  common  card  the  program  can 
be  changed  to  use  more  or  less  storage. 

/BINS/  - Bin  structure  information. 

RBIN(20)  - scoring  radii  in  (centimeters),  NRBIN  values. 

EBIN(30)  - energy  bin  boundaries  (rest  mass  units),  NEBIN+1 
values. 

UBIN(IO)  - cosine  bin  boundaries,  NUBIN  values. 

TBIN(20)  - time  bin  boundaries. 

UNCOLD(20)-  uncollided  number  flux  contribution  at  each  scoring 
radii. 

DELE(30)  - 1/absolute  value  of  the  energy  bin  width. 

DELT(20)  - 1/absolute  value  of  the  time  bin  width. 

DELU(IO)  - 1/absolute  value  of  the  angle  bin  width. 


/CDATA/  - Input  variables. 

RCALC(20)  - input  radii  converted  to  centimeters. 

NIOP  - number  of  folding  operations  requested. 

NR  - number  of  RCALC  values. 

RINPUT(20)-  input  radii  converted  to  centimeters. 

PIRSQ  - = 0 if  values  are  not  multiplied  by  4ttR^ 

= 1 if  values  are  multiplied  by  4?rR  . 

ARHOL(7,  20)  - Hollerith  identifier  associated  with  user  supplied 
angular  response  functions  ARESP1  to  ARESP7. 

ERHOL(7, 20)  - Hollerith  identifier  associated  with  user  supplied 
energy  response  functions  ERESP1  to  ERESP6. 

Y - source  yield,  default  value  1.0  photons. 

ITSOR(25)  - = 1 if  time  convolution  desired 

= 0 if  not. 

TSORT(50)  - user  supplied  time  values  for  time  convolution. 

TSORW(50)  - user  supplied  weight  associated  with  TSORT 
for  time  convolution. 

NTSOR  - number  of  TSORT, TSORW  pairs. 

/CESOR/  - Information  about  the  source  energy  spectrum. 

ESOR  - source  energy  normalization. 

ESBAR  - average  energy  of  the  source  spectrum. 

SORN(27)  - normalized  source  energy  spectrum. 

ESORKT  - ESOR,  converted  to  KT  of  photons. 

YKT  - Y,  converted  to  KT  of  photons. 

IXORG  - = 1 for  x-ray  input  tape 

= 2 for  gamma  ray  input  tape. 
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/DIGIT/  - Hollerith  constants  for  input  handling  subroutines. 

IDG(IO)  - Hollerith  representation  of  integers  0-9. 

/FOLDER/  - Information  needed  to  fold  and  output  parameters  requested. 

IFOLD(25)  - index  defining  each  type  of  fold  operation  re- 
quested. 

1 = no  folding  requested 

2 = fold  on  angle  parameter 

3 = fold  on  time  parameter 

4 = fold  on  energy  parameter 

5 = fold  on  angle  and  time 

6 = fold  on  angle  and  energy 

7 = fold  on  energy  and  time 

8 = fold  on  angle,  time  and  energy. 

IFOLDS(25)  - the  number  of  storage  locations  in  BC  required 

by  each  fold  operation. 

LFOLD(25)  - location  of  first  value  of  each  operation  in  BC 

relative  to  the  starting  location  of  each  output 
radii  requested. 

IFLAG(25)  - index  indicating  the  type  of  output  requested  for 

each  operation. 

1 = print  delta  data 

2 = print  group  data 

3 = print  normalized  delta  data 

4 = punch  delta  data 

5 = punch  group  data 

6 = punch  normalized  delta  data 

7 = write  delta  data  tape 

8 = write  group  data  tape 

9 = write  normalized  delta  data  on  tape. 

IR  - index  of  source  radii  associated  with  a particular 

flux  differential. 

IE  - index  of  energy  bin  associated  with  a particular 

flux  differential. 
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IT 

index  of  time  bin  associated  with  a particular 
flux  differential. 

IMU 

- 

index  of  cosine  bin  associated  with  a particular 
flux  differential. 

IANG(25) 

- 

index  of  angle  response  function  requested  for 
each  operation;  1-3  refer  to  built  in  response 
functions,  5-10  refer  to  input  supplied  functions. 

NTNMU 

- 

NTBIN*NUBIN  (number  of  time  bins  * number 
of  cosine  bins). 

NTNE 

- 

NTBIN*NEBIN. 

NTNUNE 

- 

NUBIN*NEBIN. 

/NGXI/ 

ICC 

- 

used  by  ATR  subroutines. 

NGX 

- 

used  by  ATR  subroutines. 

ITCH 

- 

pointer  to  which  character  in  ICH  is  currently 
being  processed. 

ICH(80) 

- 

80  column  input  card  image. 

/OUTP/ 

Input  and  output  unit  numbers. 

NONPG 

“ 

specifies  maximum  number  of  output  values  on 
a print  line.  Set  at  8 by  the  program,  can  be 
set  smaller  to  produce  a narrower  output  page. 

NOUT 

- 

output  unit,  set  equal  6. 

NIN 

- 

input  unit,  set  equal  5. 
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/POINT/ 


Storage  information  for  FOLD. 


LPOINT(19) 

the  nth  value  points  to  location  in  KPOINT  as- 
sociated with  the  nth  RBIN. 

KPOENT(40) 

locations  in  BCiess  and  BCcreater  to  be  used 
by  FOLD. 

LS(18) 

the  number  of  locations  in  KPOINT  associated 
with  each  RBIN. 

PL(20) 

used  to  interpolate  final  sum  to  the  requested 
RCALC. 

PG(20) 

used  to  interpolate  final  sum  to  the  requested 
RCALC. 

Further  explanation  of  LPOINT,  KPOINT,  and  LS  is  required  since 
they  direct  the  storage  of  data  in  the  blank  common  array  BC.  When- 
ever data  from  a new  radius  (IR)  is  read,  LS(IR)  is  examined  to 
determine  how  many  RCALC 's  will  need  this  data  for  interpolation. 

If  LS(IR)  is  zero,  the  data  at  this  radius  is  skipped.  LPOINT(IR)  then 
points  to  a location  in  KPOINT  where  locations  in  BC  are  stored. 

LS(IR)  locations  in  BC  are  stored  in  KPOINT  beginning  at  LPOINT(IR). 
Thus  the  data  will  be  stored  in  BC  LS(IR)  times,  and  the  locations  of 
these  stores  are  stored  in  the  KPOINT  array  beginning  at  LPOINT(IR). 

/RESP/  - Response  functions. 

SOR(l,  27) 

input  on  *ESOR  card. 

RANG(10,  7) 

ten  angular  response  functions  of  NUBIN+1 
values  each.  Functions  1 through  3 are  built 
into  the  code,  4-10  may  be  input. 

RENG(10,  27) 

ten  energy  response  functions  with  NEBIN  values 
each.  Functions  1 through  4 are  preset,  5-10  may 
be  input. 
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/TDATA/ 


Contains  data  read  from  input  tape. 
- source  energy  number. 


ISOR 
NRBIN 
NEBIN 
NTBIN 
NUBIN 
LENGTH 

IBUF 

NREC 

ESORHI 

ESORLO 


number  of  radii  for  source  energy  band. 

number  of  energy  bins. 

number  of  time  bins. 

number  of  angle  bins. 

number  of  data  for  this  source  band 
= NRBIN  *(NE  BIN  -LSOR+ 1)  *NTBIN  *NUBIN . 

length  of  FLAIRIT  input  record  = 5000. 

number  of  records  necessary  to  transmit  LENGTH 
words.  The  last  record  may  be  a partial  record. 

upper  energy  of  source  band  (rest  mass  units). 

lower  energy  of  source  band  (rest  mass  units). 


3.  2 FLAIR  SUBROUTINE  DESCRIPTION 
Subroutine  CARDIN 

This  subroutine  initializes  and  selects  default  values  for  required 
program  variables.  A test  is  made  to  determine  the  type  of  card 
being  read,  then  a branch  executed  to  the  appropriate  code  for  pro- 
cessing. 


The  *R  card  is  checked  for  units.  If  none  have  been  specified,  centi- 
meters are  assumed.  The  values  are  read  and  converted  to  centi- 
meters if  necessary.  More  than  one  *R  card  may  be  used  to  input 
a maximum  of  20  radii. 


The  *0  card  is  searched  for  various  parameters  which  may  be  in 
any  order.  Default  values  are  set  for  any  that  are  missing.  A 
blank  *0  card  will  be  interpreted  as  *0  PH  CUR  TIME  D.  Any 
duplicate  parameters  on  the  card  will  be  ignored.  The  *T  and  *P 
cards  are  processed  in  the  same  manner  as  a *0  card,  except  tape 
or  punch  output  is  flagged. 

The  *ARESPn  card  is  searched  for  a title  (maximum  20  characters) 
flagged  with  a * before  and  after.  If  none  is  found  ARESPn  will  be 
used.  NUBIN  or  NUBIN+1  values  are  required.  The  NUBIN+1  value 
is  applied  to  the  uncollided  component.  If  it  is  omitted  it  will  be  set 
equal  to  the  first  value.  If  less  than  NUBIN  values  are  read,  the 
card  will  be  flagged  as  illegal.  A maximum  of  7 *ARESP  cards 
may  be  used. 

The  *ERESPn  card  is  also  searched  for  a title  (maximum  20  char- 
acters). If  no  title  has  been  punched,  ERESPn  will  be  used.  The 
program  expects  NEBIN  values.  If  too  few  values  are  found,  the 
card  will  be  considered  illegal.  A maximum  of  6 *ERESP  cards 
may  be  used. 

The  *ESOR  card  is  checked  for  units.  If  no  units  are  specified 
photons/group  are  assumed.  The  values  are  read  and  converted 
to  program  units  if  necessary.  If  less  than  NEBIN  values  are  input, 
the  array  will  be  zero  filled.  These  values  are  summed  to  obtain 
the  source  normalization, and  then  used  to  compute  average  energy  of 
the  source  spectrum.  Only  one  *ESOR  card  is  used  per  *EXC  card. 
Duplicate  *ESOR  cards  will  be  ignored. 

The  *DEN  card  is  checked  for  units.  If  the  units  have  been  left 
blank,  GM/CM3  are  assumed.  The  value  read  from  the  card  is 
converted  to  program  units  if  necessary.  If  no  value  is  found,  the 
card  will  be  ignored.  If  the  *DEN  card  is  not  encountered  the  code 
uses  the  default  value  of  1.  29  x 10-3  GM/CM3.  If  more  than  one 
♦DEN  is  found  it  will  be  ignored. 


The  *TSOR  card  is  read  and  checked  for  an  even  number  of  values. 
The  program  expects  time-weight  pairs,  i.e. , 

*TSOR  Time  1 Weight  1 Time  2 Weight  2 ... 

Any  additional  *TSOR  cards  are  ignored. 
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The  *Y  card  is  checked  for  unit  specification.  If  it  has  been  omitted, 
units  of  photons  are  assumed.  One  value  is  expected  by  the  code.  If 
one  is  not  found,  the  card  will  be  considered  illegal  and  be  ignored. 
The  value  that  is  found,  is  converted  to  program  units  and  is  used 
instead  of  the  source  normalization  computed  from  the  *ESOR  card. 
The  *Y  card  is  optional  and  only  one  should  occur  before  each  *EXC 
card.  Additional  ones  will  be  ignored. 

The  *EXC  signals  end  of  all  cards  needed  to  process  one  batch  job. 

If  4PIRSQ  appears  on  the  card,  4trR2  values  will  be  output. 

The  **  card  signals  end  of  job.  The  size  of  BC  limits  the  number  of 
operations  that  can  be  processed  at  one  pass  as  a function  of  number 
of  radii  and  type  of  operation.  Multiple  batches  can  be  processed  on 
one  run.  If  the  **  card  does  not  follow  the  *EXC  card  the  code  at- 
tempts to  process  another  batch. 

After  the  *EXC  has  been  read,  CARDIN  checks  to  see  if  enough  in- 
formation has  been  input  to  proceed.  If  required  information  is 
missing  and  no  default  can  be  supplied,  the  program  stops  with  an 
error  message. 

A synopsis  of  the  folding  operations  requested  is  printed  to  tell  the 
user  how  the  code  has  interpreted  the  input  cards.  If  all  criteria  to 
proceed  have  been  satisfied  control  returns  to  the  main  program. 


Subroutines  GETVAL,  RDWRTC,  FREFNO,  IN  TERR,  INGET, 
NNBCHR  were  borrowed  from  ATR  to  be  used  by  CARDIN.  These 
routines  convert  Hollerith  data  to  floating  point  or  integer  values. 

Called  from:  FLAIR. 


J . I , $ f ' 


Subroutines  called:  GETVAL,  RDWRTC,  NNBCHR,  NCHR. 

Common*  required:  Blank,  ACCSRY,  BIGGER,  BINS,  CDATA, 

CESOR,  DIGIT,  FOLDER,  IRDINK,  NGXI, 
OUTP,  RESP. 


Variables  required:  See  input  instructions  in  Appendix  A and  through- 
out the  text  of  Section  2. 
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Main  Program  FLAIR 


! 


FLAIR  sets  the  storage  limit  for  BC  and  initializes  it  to  zero.  The 
first  record  of  the  data  tape  (tape  3)  is  read  to  determine  if  it  is  an 
x-ray  tape  or  gamma -ray  tape. 

The  first  record  or  header  record  associated  with  ISOR-1  is  read 
to  determine  bin  structure.  The  header  record  for  succeeding  source 
groups  will  be  ignored  except  for  values  which  change,  i.e. , ISOR, 
LENGTH,  NREC,  ESORHI,  ESORLO,  RBIN  and  UNCOLD). 

Subroutine  CARDIN  is  called  to  read  and  process  all  input  cards  up 
to  and  including  a *EXC  card.  An  error  in  input  will  cause  a pro- 
gram stop  in  CARDIN. 

If  source  energy  group  ISOR  has  been  requested  by  the  *ESOR  card, 
RTEST  is  called  to  determine  which  source  energy  radii  are  needed 
to  compute  values  of  BC^ess  and  BC,rreater.  The  data  corresponding 
to  each  flagged  RBIN  is  processed,  other  lata  is  ignored.  A call  to 
UNPR  gives  values  of  IE,  IT,  IMU,  IR  for  each  data  value.  FOLD 
adds  the  value  to  appropriate  location  in  BCiess  or  BCgrea^er;  When 
all  records  associated  with  source  group  ISOR  have  been  read,  the 
uncollided  component  is  added  through  a call  to  entry  UNCOLL. 

All  values  in  the  BCiess  and  BCgreater  arrays  are  modified  by  the 
source  normalization  or  by  the  value  input  on  the  *Y  card.  SUM  is 
called  to  interpolate  between  BCiess  and  BC  greater  and  add  to  the 
accumulative  totals  in  BCsum.  Each  required  source  group  is 
handled  in  the  same  manner  until  all  required  groups  are  processed. 

If  4ttR2  values  have  not  been  requested  , BCSUm  is  divided  by  4wr2. 

If  time  convolution  is  requested  TIMCON  is  called.  If  x-ray  data  is 
being  input  EBIN  is  converted  to  keV  for  output.  After  a call  to 
OUTPUT,  the  input  tape  is  rewound  and  control  returns  to  the  begin- 
ning to  work  on  the  next  batch  if  one  has  been  requested. 

Subroutines  called:  CARDIN,  RTEST,  UNPR,  FOLD,  SUM, 

TIMCON,  OUTPUT. 
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Commons  required:  Blank,  BINS,  FOLDER,  RESP,  TDATA,  CDATA, 

POINT,  OUTP,  BIGGER,  CESOR. 

Units  required:  INPUT,  OUTPUT,  TAPE5=INPUT,  TAPE6= 

OUTPUT,  TAPE  3. 

Optional  units:  TAPE 7,  TAPE8. 

Subroutine  FOLD  (IOP,  IWAY,  DATA) 

Subroutine  FOLD  processes  one  value  of  differential  flux  (DATA)  on 
each  call.  Knowing  the  corresponding  angle  bin  (IMU)  ar.d  energy 
bin  (IE),  DATA  is  multiplied  by  the  proper  angle  and  energy  response 
function.  Using  the  pointers  determined  by  RTEST  and  knowing  the 
type  of  folding  operation  requested,  the  proper  storage  locations  in 
BCiess  or  BCgreater  can  be  determined.  DATA  is  then  added  to  the 
computed  locations. 

The  entry  UNCOLL  is  used  to  add  the  uncollided  component  to  the 
proper  location,  after  being  modified  by  the  angle  and  energy  response 
functions.  See  the  layout  of  BCiess  and  BCgreater  in  the  blank  common 
description. 

Called  from:  FLAIR. 

Subroutines  called:  none. 

Entry  points:  UNCOLL. 

Commons  required:  Blank,  FOLDER,  POINT,  RESP,  TDATA. 
Variables  required: 

IOP  - index  of  current  folding  operation. 

IWAY  - index  for  a computed  GO  TO 

1 = total  operation  (no  folding) 

2 = angle  dependent  operation 

3 = time  dependent  operation 

4 = energy  dependent  operation 
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angle  and  time  dependent  operation 
angle  and  energy  dependent  operation 
time  and  energy  dependent  operation 
angle,  time  and  energy  dependent  operation 
total  operation  uncollided 
angle  and  time  dependent,  time  convoluted, 
uncollided 

time  dependent,  time  convoluted,  uncollided 
time  and  energy  dependent,  time  convoluted, 
uncollided 

angle,  time  and  energy  dependent,  time  con- 
voluted, uncollided 

DATA  - the  current  four-dimensional  differential  flux  value. 

Variables  changed:  each  call  to  FOLD  changes  all  locations  in 

Euless  and  ^greater  requiring  DATA. 

Subroutine  OUTPUT  (IOP,  LE,  L) 

Subroutine  OUTPUT  offers  four  printed  output  options;  group  data, 
normalized  group  data,  delta  data  and  normalized  delta  data.  If 
an  option  has  not  been  specified  on  the  *0  card  read  by  CARDIN, 
delta  data  will  be  selected. 

If  group  data  is  uesired  the  BCsum  array  is  output  as  is.  If 
normalized  group  data  is  flagged,  each  value  in  BCsum  for  this 
fold  operation  will  be  divided  by  the  sum  of  these  values,  and  then 
output.  If  delta  data  is  requested,  each  value  in  BCsum  as- 
sociated with  the  fold  operation  is  divided  by  the  appropriate  delta 
values.  If  normalized  delta  values  are  requested,  the  appropriate 
BCSUm  array  is  divided  by  the  array  sum  and  by  the  required  delta 
values. 

Appropriate  titles  are  selected  to  define  the  output  to  the  user. 

One  fold  operation  is  output  with  each  call  to  OUTPUT.  The 
maximum  number  of  columns  output  per  page  is  set  by  NONPG. 

The  code  sets  N0NPG=8,  this  can  be  reduced  to  produce  a narrower 
page  if  desired.  If  an  entire  line  of  output  is  all  zeros,  it  will  not 
be  printed. 

Card  or  tape  output  are  available  through  a call  to  OUT 2.  Printed 
output  will  also  be  produced. 
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Called  from:  FLAIR. 


Subroutines  called:  OUT2. 

Commons  required:  Blank,  BINS,  CDATA,  CESOR,  DIGIT, 

FOLDER,  OUTP,  TDATA. 

Variables  required: 

IOP  - index  indicating  the  fold  operation  being  output. 

LE  - starting  location  in  BCsum  for  operation  IOP 

at  the  first  radius. 

L - the  number  of  locations  associated  with  one 

radius . 

Variables  changed:  If  other  than  group  data  is  requested,  the  storage 

locations  in  BCsum  associated  with  operation  IOP 
are  changed  as  requested. 

Subroutine  OUT2  (ITAPE,  IOP,  LE,  L,  IFOLDS) 

If  punched  output  is  requested,  the  first  card  contains  IOP,  RCALC(I), 
IFOLDS.  Format  is  15,  5X,  E10.4,  110.  Then  IFOLDS  values  from 
BCsum  are  punched  with  format  8E10.  4.  This  pattern  is  repeated 
for  each  RCALC. 

If  a tape  is  desired,  2 records  of  unformatted  data  are  written  for 
each  RCALC.  The  first  record  contains  IOP,  RCALC(I),  IFOLDS. 

The  second  record  contains  IFOLDS  values  from  BCsum  associated 
with  IOP  and  RCALC(I). 

Called  from:  OUTPUT. 

Subroutines  called:  none. 

Commons  required:  CDATA. 
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Variables  required: 


ITAPE  - file  to  be  punched  or  written. 

IOP  - index  of  current  fold  operation. 

LE  - starting  location  of  BCgum. 

L - number  of  locations  associated  with  each  RCALC. 

IFOLDS  - number  of  storage  locations  to  output. 

Subroutine  RTEST 

Subroutine  RTEST  is  called  once  for  each  source  energy  band  to  be 
processed.  The  code  determines  which  source  energy  band  radii 
(RBIN(J))  will  be  needed  for  storing  values  in  BCiess  and  BCgreater* 

It  further  decides  how  many  times  RBIN(J)  associated  values  will  be 
stored  and  the  starting  locations  for  each  store. 

Values  for  PG(I)  and  PL(I)  are  computed  for  each  requested  input 
radius  (RCALC(I))  to  be  used  later  by  SUM  to  interpolate  the  results 
between  BC^ess  and 

If  RCALC(I)  is  less  than  RBIN(l),  nothing  will  be  stored  in  BC^ss* 
Later  the  BCgreater  value  will  be  added  to  BCsum.  If  RCALC(I) 
is  greater  than  RBIN(NEBIN),  nothing  is  stored  in  BCjess  or 
BCgreater*  Later  a value  of  zero  will  be  added  to  BCsum. 

Called  from:  FLAIR. 

Subroutines  called:  none. 

Commons  required:  BINS,  CDATA,  FOLDER,  POINT,  TDATA. 
Variables  required:  none. 

Variables  changed: 

LPOINT(19)-  LPOINT(J)+l  points  to  the  first  location  in 

KPOENT  to  find  the  first  BC  storage  location 
for  data  associated  with  RBIN(J). 


L 
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KPOENT(40)  - 

starting  locations  in  BCiess  and  BC^eater 
to  be  used  by  FOLD. 

LS(18) 

LS(J)  gives  the  number  of  locations  in  KPOINT 
associated  with  RBIN(J). 

PL(20) 

power  applied  to  BCiess  during  interpolation 
process.  One  value  for  each  RCALC. 

PG(20) 

power  applied  to  BCgreater  during  interpolation 
process.  One  value  for  each  RCALC. 

Subroutine  SUM 

After  all  required  data  has  been  read  for  one  source  energy  group, 
SUM  is  called  to  interpolate  between  BCiess  and  BCgreater-  If 
PL  equals  zero,  no  interpolation  is  required,  and  the  BCg^a^ 
value  is  added  to  BCsum.  If  PL  is  non  zero,  a semilog  interpolation 
is  used  to  determine  the  value  added  to  BCsum. 

Called  from:  FLAIR. 

Subroutines  called:  none. 

Commons  required:  Blank,  CDATA,  FOLDER,  POINT. 

*■-  Variables  required:  none. 

Variables  changed:  all  values  of  BCsum;  BCieSs  and  BC^pa^-  are 

set  to  zero  upon  return  to  FLAIR. 

Subroutine  TIMCON 

Subroutine  TIMCON  is  called  if  time  convolution  has  been  requested. 
The  values  of  TSORT  and  TSORW  are  processed  to  create  a matrix  M. 
The  type  of  fold  operation  is  determined  in  order  to  properly  set  the 
DO  index  and  incrementing  parameters  needed  to  access  correct  values 
in  BCsum-  One  fold  operation  is  handled  at  a time.  All  collided  com- 
ponents are  multiplied  by  matrix  M and  summed  in  vector  C.  The 
BCsum  values  are  then  replaced  by  C.  The  uncollided  component  is 
added  to  the  proper  location  in  BCSum- 
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Called  from:  FLAIR. 


Subroutines  called:  none. 

Commons  required:  Blank,  BINS,  CDATA,  FOLDER,  TDATA. 
Variables  required: 

J1  - Jl+1  is  location  of  first  word  of  BCsum. 

L - number  of  locations  associated  with  one  RCALC. 

Variables  changed:  locations  in  BCsum  which  should  be  time 

convoluted. 

Subroutine  UNPR(I) 

Subroutine  UNPR  computes  or  "unpacks"  the  indices  IR,  IE,  IT,  and 
IMU  for  the  Ith  differential  flux.  I was  computed  in  DART  using  the 
following  equation 

I = IE  + NE  BIN  (IT-1+NTBIN*(IMU-1+NUBIN*(IR-1))) 

where  IE  runs  from  ISOR  to  NEBEN. 

The  subroutine  checks  if  the  correct  value  of  I is  one  greater  than 
the  last  one  processed.  If  it  is,  a short  cut  is  used  to  determine 
IR,  IE,  IT  and  IMU. 

Called  from:  FLAIR. 

Subroutines  called:  none. 

Commons  required:  FOLDER,  TDATA. 

Variables  required: 

I - index  to  be  unpacked,  varies  from  1 to  LENGTH. 
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Variables  changed: 


IR  - index  of  associated  radii. 

IE  - index  of  associated  energy  bin. 

IT  - index  of  associated  time  bin 

IMU  - index  of  associated  cosine  bin. 
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APPENDIX  A 


SUMMARY  OF  FLAIR  INPUT  INSTRUCTIONS 


The  following  card  types  are  available  in  FLAIR: 

♦R 

*0 

*T 

*P 

♦DEN 

♦ESOR 

♦Y 

♦TSOR 

♦ARESPn 

♦ERESPn 

♦EXC 

** 


The  ♦ is  punched  in  column  1,  the  identifier  follows  immediately 
with  no  blanks  allowed. 

Each  run  requires  at  least  one  *R  card,  at  least  one  *0  card  (or 
♦T  or  *P),  an  *ESOR  card,  and  an  *EXC  card. 

The  parameters  on  each  card  are  separated  by  one  or  more  blanks. 
Parameters  may  be  continued  on  a second  card,  except  for  the  *0,  *T, 

*P  cards. 

Range  Card  - Sets  the  ranges  where  results  will  be  calculated. 

*R  units  values 

units  may  be  KM,  KFT,  MI,  FT,  M,  CM  (default),  GM 

More  than  one  *R  card  may  be  used.  Maximum  number  of  radii 
per  rim  is  20. 
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Density  Card  - Sets  the  air  density. 
*DEN  units  values 


units  may  be  GM/CM3  (default),  AT/BRNCM 

If  *DEN  card  is  omitted  a value  of  1.  29  x 10  "3  gm/cm3  will  be 
used. 

Source  Spectrum  Card  - Defines  the  energy  dependence  of  the  source. 
*ESOR  units  values 

units  may  be  PH  (default),  MEV,  KEV,  CAL,  PH/ME V,  PH/KEV, 
E/E,  CAL/KEV 

One  *ESOR  card  is  required  for  each  execution.  Any  normalization 
that  exists  on  the  *ESOR  card  is  over-ridden  by  a yield  card. 

Yield  Card  - Sets  the  source  normalization. 

*Y  units  values 

units  may  be  PH  (default),  KT,  MEV,  CAL,  MT,  MO 

If  used,  a yield  card  over-rides  the  source  normalization  used  on 
the  *ESOR  card. 

Output  Card  - Specifies  the  output  including  responses,  dimensionality  of 
results  and  form  of  output. 


*0  angle  energy  time  type  flag 

These  parameters  may  occur  in  any  order. 

angle  may  be  or  angle  dependent 


FLUX 

FCUR 

CUR  (default) 

ARESPn  (7  is  max.  value  of  n) 


DFLUX 
DFCUR 
DC  UR 
DARES  Pn 
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energy  may  be 


or  energy  dependent 


PH  (default)  DPH 

MEV  DMEV 

KEV  DKEV 

CAL  DCAL 

ERESPn  (6  is  max.  value  of  n)  DERESPn 


time  may  be 


or  time  dependent 


TIME  (default) 


DTI  ME 
DTSOR 


type  may  be 

G for  group  data  output 
GN  for  normalized  group  data  output 
D for  delta  data  output  (default) 

DN  for  normalized  delta  data 


flag  may  be  blank  or  ONLY. 

If  ONLY  is  specified,  only  the  operation  requested  will  be  output. 

If  left  blank, all  lower  level  operations  will  be  output  as  well.  For 
example,  if  a time  and  angle  dependent  operation  is  requested,  the 
time  dependent  output,  the  angle  dependent  output,  and  the  total 
data  will  be  output. 

Tape  Card  - Specifies  output  to  tape. 

*T  angle  energy  time  type  flag 

This  is  the  same  as  an  *0  card,  but  data  will  be  written  on  tape 
unit  7. 

Punch  Card  - Specifies  output  to  punch. 

*P  angle  energy  time  type  flag 

This  is  the  same  as  an  *0  card,  but  data  will  be  written  on  unit  8 
for  punching. 
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Angle  Response  Card  - Sets  user  input  angle  response. 

* ARESPn  *title*  values 
n = 1,  7 (maximum  of  7 angle  response  cards) 

title  - is  a maximum  of  20  characters  enclosed  by  asterisks, 
if  omitted  ARESPn  will  be  used. 

values  - one  value  needed  for  each  angle  bin,  plus  one  more 
for  uncollided  data.  If  last  value  not  input  it  will  be 
set  equal  to  the  first  value. 

This  card  required  when  ARESPn  or  DARESPn  is  specified  on  *0, 

*T  or  *P  cards. 

Energy  Response  Card  - Sets  user  input  energy  response. 

*ERESPn  *title*  values 

n = 1, 6 (maximum  of  6 energy  response  cards) 

title  - is  a maximum  of  20  characters  enclosed  by  asterisks, 
if  omitted  ERESPn  will  be  used. 

values  - one  value  needed  for  each  energy  bin. 

This  card  required  when  ERESPn  or  DERESPn  is  specified  on  *0, 

*T  or  *P  cards. 

Execution  Cards  - Indicates  completion  of  input. 

*EXC  or  *EXC  4PIRSQ 

The  *EXC  card  signals  the  end  of  input  batch.  If  4PIRSQ  is  specified 
all  values  are  multiplied  by  4ffR2  where  R is  the  radius  in  centimeters. 

End  of  Job  Card  - Indicates  end  of  job. 

** 

The  **  card  signals  end  of  job,  allowing  multiple  batches  each  ending 
with  the  *EXC  card.  If  **  card  is  missing  job  will  terminate  trying 
to  read  the  next  batch  after  finishing  the  previous  batch. 
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APPENDIX  B 

THREE  SAMPLE  FLAIR  CALCULATIONS 
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This  Appendix  presents  three  sample  FLAIR  calculations  that  were 
selected  to  demonstrate  both  the  use  and  utility  of  the  FLAIR  code.  Tne 
first  two  are  x-ray  calculations  with  a 15  kT  blackbody  spectrum  and  the 
third  is  a gamma -ray  calculation.  To  execute  FLAIR  either  the  x-ray  or 
the  gamma -ray  data  tape  which  ever  is  appropriate  must  be  mounted  on 
unit  three.  The  input  cards  used  for  the  sample  problems  are  listed  as 
part  of  the  problem  output.  Each  of  the  three  problems  executes  in  less 
than  14  seconds  on  a CDC  7600  computer  for  a cost  of  less  than  five  dollars 
each. 

The  first  calculation  computes  forward  current  calories  from  a 15  kT 
blackbody  source  at  thirteen  ranges,  giving  both  total  calories  and  energy 
spectra.  The  second  sample  calculation  demonstrates  time  convolution  as 
well  as  other  options.  Again  the  x-ray  source  is  a 15  kT  blackbody  this  time 
at  a 1000  ft  range.  Time  dependence  is  calculated  with  all  the  source  at 

Q 

time  zero  and  with  the  source  spread  over  two  shakes  (2  x 10 ~ seconds). 
Energy  spectra  and  angle  dependence  are  also  calculated.  The  third  cal- 
culation is  a gamma -ray  calculation  that  demonstrates  the  use  of  the  user 
input  energy  and  angle  response  functions.  The  assumed  response  functions 
represent  dose  in  rads(Si)  in  a missile  vulnerability  calculation,  although 
the  response  function  themselves  are  ficticious  and  do  not  relate  to  any 
existing  missile.  The  gamma -ray  source  spectrum  also  is  ficticious  and 
is  to  be  used  only  as  a sample  problem.  Fluences  and  sample  missile  re- 
sponses are  calculated  at  three  ranges. 
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No.  of  No.  of 

Copies  Organization  Copies  Organization 


12  Commander 

Defense  Documentation  Center 
ATTN:  DDC-TCA 

Cameron  Station 
Alexandria,  VA  22314 

1 Director 

Defense  Nuclear  Agency 
ATTN : RATN 

Washington,  DC  20305 

1 Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCDMA-ST 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

1 Commander 

US  Army  Aviation  Systems 
Command 

ATTN:  DRSAV-E 

12th  and  Spruce  Streets 

St.  Louis,  MO  63166 

1 Director 

US  Army  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

1 Commander 

US  Army  Electronics  Command 

ATTN:  DRSEL-RD 

Fort  Monmouth,  NJ  07703 

1 Commander 

US  Army  Missile  Command 

ATTN:  DRSMI-R 

Redstone  Arsenal,  AL  35809 


1 Commander 

US  Army  Tank  Automotive 
Development  Command 
ATTN:  DRDTA-RWL 

Warren,  MI  48090 

2 Commander 

US  Army  Mobility  Equipment 

Research  § Development  Command 
ATTN:  Tech  Docu  Cen , Bldg.  315 

DRSME-R2T 

Fort  Belvoir,  VA  22060 
1 Commander 

US  Army  Armament  Command 
Rock  Island,  IL  61202 

1 Commander 

US  Army  Frankford  Arsenal 
ATTN:  SARFA-FCV,  Stan  Goodman 

Philadelphia,  PA  19137 

2 Commander 

US  Army  Picatinny  Arsenal 
ATTN:  SARPA-ND-C-T 

JAWTIP/SARPA-RT-S 
Dover,  NJ  07801 

1 Commander 

US  Army  Harry  Diamond  Labs 
ATTN:  DRXDO-TI 

2800  Powder  Mill  Road 
Adelphi , MD  20785 

1 Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SA 

White  Sands  Missile  Range 
NM  88002 
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1 Commander 

US  Army  Nuclear  Agency 
Fort  Bliss,  TX  79916 

1 AFWL  (SUL) 

Kirtland  AFB,  NM  87117 

1 Director 

Lawrence  Livermore  Laboratory 
ATTN:  Tech  Info  Div 

P.  0.  Box  808 
Livermore,  CA  94550 

1 Director 

Los  Alamos  Scientific  Lab 
ATTN : Rept  Lib 

P.  0.  Box  1663 
Los  Alamos,  NM  87544 

2 Science  Applications,  Inc. 

ATTN:  W.  H.  Scott,  Jr. 

B.  L.  Colborn 
1200  Prospect  Street 
La  Jolla,  CA  92037 

Aberdeen  Proving  Ground 

Marine  Corps  Ln  Ofc 
Dir,  USAMSAA 
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